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We studied the effects of the incommensurate structural modulations on the ladder subsystem of 
the Sr\4,- x Ca x Cu2iO a family of compounds using ab-initio explicitly-correlated calculations. From 
these calculations we derived t — J model as a function of the fourth crystallographic coordinate r 
describing the incommensurate modulations. It was found that in the highly calcium-doped system, 
the on-site orbital energies are strongly modulated along the ladder legs. On the contrary the two 
sites of the ladder rungs are iso-energetic and the holes are thus expected to be delocalized on the 
rungs. Chain-ladder interactions were also evaluated and found to be very negligible. The ladder 
superconductivity model for these systems is discussed in the light of the present results. 

PACS numbers: 71.10.Fd, 71.27. +a, 71.23.Ft 
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I. INTRODUCTION 

The Sru- x Ca x Ou2iOn family of compounds has at- 
tracted a lot of attention in the last years due to the 
discovery of a superconducting phase 1 in these quasi- 
unidimcnsional systems. The Sri4\- x Ca x Cu2iO^\ are lay- 
ered materials with alternated planes of two-leg spin lad- 
ders and planes of spin chains 2 . The superconducting 
phase, observed in the highly calcium-doped systems, at 
low temperature and under pressure, is believed to be the 
realization of the predicted superconductivity in two-leg 
ladders systems^. Indeed, it is currently assumed that the 
Sri i ^ x Ca x Cki2iOu ladder subsystem can be represented 
by a two-leg ladder isotropic t — J model with J/t ~ 0.4 
and a hole doping ranging between 1— and 2.8^ holes per 
formula unit (f.u.), that is between 0.07 and 0.2 hole per 
ladder site. Numerical as well as slave-bosons analyti- 
cal calculations 6 , showed that in this parameter range, 
the superconducting pairing fluctuations are dominant. 
Hole-pairs are predicted to be bounded on the rungs and 
the superconductivity to be supported by pairs collective 
modes. 

The Sri^- x Ca x Ou,2iOn compounds present structural 
incommensurate modulations of both the chains and lad- 
der subsystems. It was showni using ab-initio calcula- 
tions, that, unlike what was usually assumed, the chain 
subsystem cannot be considered as a quasi-homogeneous 
system for which the structural modulations induce only 
small perturbative effects. In fact, the modulations 
caused by the influence of the ladders on the chain sub- 
system are responsible for the low-energy physics of the 
latter, that is charge localization and spin arrangement. 
These ab-initio calculations showed that the parameters 
of the second neighbor t — J + V model are strongly mod- 



ulated. The major effects are (i) on the on-site magnetic 
orbital energies that fluctuate within a few electron- Volt 
range, (ii) on the first neighbor effective exchange in- 
tegrals that range from ferromagnetic (as expected for 
such 90° oxygen-bridged copper sites) to antiferromag- 
netic (in the calcium highly-doped compounds). More- 
over the calcium iso-electronic substitution, that acts as 
a chemical pressure, is responsible for large variations 
of the t — J model parameters. It results a quite differ- 
ent chain ground-state for the undoped and highly-doped 
systems. The undoped compound exhibits the formation 
of second-neighbor dimers caused by the electron local- 
ization on sites with low orbital energy, while the large 
x compounds exhibit the formation of low spin clusters 
with antiferromagnetic local ordering. 

One can thus wonder whether the structural modula- 
tions induce as important effects on the ladder electronic 
structure, as observed on the chain one. This point is 
of crucial importance since the existence of such an ef- 
fect would strongly questioned the interpretation of the 
observed superconductivity. Another point, often ad- 
dressed in the literature, is the importance of possible 
chain-ladder coupling. Indeed, in the highly calcium- 
doped compound, some of the chain oxygen atoms come 
close to the ladder coppers as apical atoms. The aim 
of the present work is to provide answers to these ques- 
tions, by the means of ab-initio calculations. The method 
used is an embedded fragment spectroscopy method that 
explicitly computes both the correlation effects and the 
screening effects. A t — J model, explicitly dependent 
on the structural modulations, is derived from the ab- 
initio results for the ladders. The calculations are done 
on both ladders and mixed, chain and ladder, fragments, 
in the low-temperature phases of both the undoped^ and 
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x = 13.6 calcium-doped 9 compounds. It would have been 
of great interest to also run calculations in the super- 
conducting phase, unfortunately structural data are not 
available under pressure. Nevertheless we will discuss 
possible extrapolation of the ambient pressure results. 

The next section will be devoted to the computational 
details. Section 3 will discuss the t — J model, section 4 
will discuss the chain-ladder interactions and finally the 
last section will be devoted to discussions and conclu- 
sions. 



II. METHOD AND COMPUTATIONAL AND 
DETAILS 

It is well known that, in strongly correlated systems, 
few electrons (per unit cell) are responsible for their spec- 
tacular properties. These electrons are usually unpaired 
and localized both spatially and energetically near the 
Fermi level (from now on we will refer to them as mag- 
netic electrons). 

Consequently, the electronic wave function of these sys- 
tems is essentially multi-configurational, i.e. it cannot be 
described, even in a qualitative way, by a single Slater 
determinant. The consequence is that single-reference 
based methods such as Hartree-Fock plus perturbation 
theory or even Density Functional Theory (DFT) fail 
to correctly treat some characteristics of the electronic 
structure of such systems^. While DFT methods yield 
reasonable Fermi surfaces such as in the doped super- 
conductor copper oxides, the LDA method fails to pre- 
dict the correct insulating versus metallic behavior in a 
lot of correlated systems or effective local interactions 
such as the magnetic exchange constant J. This is for 
instance the case for the same copper oxides for which 
the Fermi surfaces were correctly predicted. Indeed, the 
LDA estimate of the magnetic exchange in the La^OaO^ 
compound ranges for instance from 600 meV to 800 meV 
according to the functional used while the experimen- 
tal measurements^ yield 135 ± 6 meV and the embed- 
ded cluster method used in the present work 145 meV— . 
LDA+U 13 or constrained LDA 14 (CLDA) methods that 
treat the correlation within the magnetic orbitals in a 
mean field scheme yield better results. For instance the 
LDA+U method yields between 129 and 165 meV for 
the J value of Sr^CuO^CL-i (exp. 125 meV) and the 
CLDA method yields 128 meV for the La-JJuO^ com- 
pound. These methods however still present a lot of ar- 
bitrariness, for instance the uncertainty on the results for 
the CLDA are of 50% according to the authors of refw 

A few density functional calculations are present in the 
literature for the present system. Arai and coworkers 15 
performed a LDA calculation on the average structure 
(without the subsystem modulations) of the Sr\±&±2AO^\ 
in quasi periodic approximation. Their results exhibit a 
metallic behavior for both the chains and ladders subsys- 
tems whereas they are experimentally found to be insu- 
lating. More recently, unpublished results from Schuster 



and Schwingenschlgl 16 found the same results when in- 
cluding the structural modulations. Let us point out that 
the LDA+U method should yield better results on such 
strongly correlated systems. Unfortunately calculations 
including the structural modulation necessitate a much 
too large quasi-periodic unit cell for the calculation to be 
tractable. 

For such problems it is therefore necessary to rely 
on wave-function, multi-reference, ab-initio, spectroscopy 
methods such as the Difference Dedicated Configurations 
Interaction 1 ^^ (DDCI) on embedded clusters. This 
method proved its ability to quantitatively reproduce 
the local interactions of a large family of magnetic sys- 
tems such as high T c copper oxides^, vanadium oxides^, 
nickel and cuprate fluorides^ 9 -, etc. Such multi-reference 
wave function methods are very powerful to treat strong 
correlation problems, however they necessitate a formally 
finite system. The use of such methods for the treatment 
of strongly correlated materials is possible since the in- 
teractions between the magnetic electrons are essentially 
local. In fact it was proved, both by quasi-degenerate 
perturbation analysis^ and numerically by increase of 
the cluster svz^^^S^. that only the bridging ligands 
and the first coordination shell of the magnetic centers 
are of importance for the magnetic, transfer or coulomb 
interactions. 

The DDCI method is a variational method specifically 
designed for the calculation of excited states involving 
open shells. It is based on the exact diagonalization of a 
selected configuration space built as follow. First a parti- 
tion of the cluster orbitals in occupied, active and virtual 
is defined. The active set of orbitals is defined to be the 
magnetic orbitals supporting the low energy physics as 
defined above. In the present work this set of orbitals is 
formed by the open 3d orbitals of the copper atoms, that 
is the 3d x 2_ z 2 orbitals for the ladder subsystem and the 
3d xz orbitals for the chain subsystem (x refers to the crys- 
tallographic a direction and z to the c direction). The 
occupied orbitals are the inner orbitals that one can ex- 
pect to be essentially doubly-occupied in the system and 
the virtual orbitals are the orbitals essentially empty in 
ground and low energy states of the system. The set of 
all determinants that can be built from all possible spin 
and charge arrangements within the active orbitals, all 
occupied orbitals being doubly-occupied and all virtual 
orbitals being empty, is called the Complete Active Space 
(CAS) and is considered as the reference space. The ex- 
act diagonalization within this space will properly treat 
both the open-shell character of the problem and the elec- 
tronic correlation within the active or magnetic orbitals. 
The simultaneous optimization of the orbitals will yield 
the best possible partition (CASSCF procedure). Once 
the reference space is defined the DDCI space to be di- 
agonalized is built from the CAS determinants and all 
single and double excitations on each of them except the 
ones that excite two occupied orbitals toward two vir- 
tual ones. Indeed, one can see easily from a perturbative 
analysis that at the second order the later configurations 
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bring only a diagonal shift and do not contribute to the 
excitation energies between the low energy states that 
essentially differ by their component on the CAS. 

The DDCI method has proved to be very reliable for 
the study of the strongly correlated compounds, since, 
the use of a CAS ensures to treat properly : 

1. the multi-determinantal and open-shell character of 
the reference wave function, 

2. the correlation between the magnetic electrons, 

while the selected doubly excited configurations allows to 
include^ 3 . 

1. the screening effects (dynamical polarization and 
correlation) on the different configurations of the 
reference wave function, 

2. the effects of the ligand to metal charge transfers 
mediating the interactions, 

3. and the screening effects on the latter configura- 
tions. 

One could be surprised by the importance we give to the 
ligand to metal charge transfer configurations, however 
their importance can be shown by a simple perturbative 
analysis of the microscopic mechanism supporting the ef- 
fective exchange and transfer integrals. This mechanisms 
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FIG. 1: Schematic representation of the microscopic mech- 
anism responsible for the through-bridge transfer (first line) 
and super-exchange (second line) effective integrals between 
the a and b magnetic orbitals. The arrows picture the micro- 
scopic hopping involved at each step between the magnetic a 
or b orbitals and the ligand one, 1. 



The effects of the rest of the crystal on the com- 
puted fragments are essentially the Madelung potential 
and the exclusion effects. Exclusion effects are treated 
using total-ion pseudo-potentials^ (TIPs) that prevent 
the quantum fragment electrons to extend in regions of 
space normally occupied by the electrons of the rest of 
the crystal. In the present work TIPs2I are used on the 
first shell of atoms surrounding the computed cluster. 
The Madelung potential is reproduced by a set of point 
charges located on the material atomic position. The 
charges were chosen so that the strontium and calcium 
atoms are di-cations (Ca 2+ and Sr 2+ ). The copper and 
oxygen charges where chosen as Ou 2+ and 0~ 2+s . The S 
correction is associated to the self-doping in holes. It was 
attributed to the oxygens since the holes were shown to 
be essentially localized on themSi. The 5 values for the 
chain and ladder subsystems were chosen in agreement 
with the experimental assumptions and our previous cal- 
culations^^, that is all the holes on the chains for the 
SK4OU24O41 system and one hole transfer on the lad- 
ders for the ^0.4^13. 6^*24041 compound. Concerning 
now the hole repartition within each subsystem, we have 
chosen an homogeneous repartition, corresponding to an 
average state in order not to induce localization effects 
due to a particular repartition. One should note that 
the main effect of a localization of the magnetic elec- 
trons will be a modification of the on-site orbitals ener- 
gies through the Coulomb repulsion. This self-consistent 
effect is taken into account in this work through the inclu- 
sion of the screened NN electron repulsion in the energy 
curves (see figure 0). Due to screening effects, larger dis- 
tance coulomb repulsion is expected to be of much weaker 
importance. 

Two types of fragments are used in the present work : 

• CU2OJ fragments for the ladder legs and the ladder 
rungs interactions (see figures^ and |2b), 

• Cu,20% for the chain-ladder interactions (see fig- 
ure^). 
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where A' and A are the ligand to metal excitation en- 
ergies (let us note that these two energies are somewhat 
different due to the presence or absence of a spectator 
electron on the metal-ligand bond), U is the Coulomb 
repulsion in the magnetic orbitals and t a i , tu are the lig- 
and to metal direct transfer integrals. One sees from fig- 
ure ^ an d equation ^ that in order to obtain an accurate 
determination of the effective transfer and exchange in- 
tegrals, one should properly evaluate the ligand to metal 
excitation energy. 





FIG. 2: Schematic representation of the computed fragments, 
a) for the ladder legs interactions, b) for the ladder rungs in- 
teractions and c) for the chain-ladder interactions. Circles 
represent the fragment atoms ; the gray circles represent oxy- 
gen atoms, while the black circles represent the copper atoms. 
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On-site orbital energies (e) and hopping (t) param- 
eters of the t — J model are extracted from the first 
doublet-doublet excitation energies, and associated wave- 
functions, of the above fragments with one magnetic 
electron. Exchange integrals (J) are deduced from the 
singlet-triplet first excitation energies on the same frag- 
ments with two magnetic electrons. A least-square fit 
method is used in order to fit the effective parameters to 
the ab-initio results. 

The calculations were done using the MOLCAS2& and 
CASDI set of programs and the basis sets used can be 
found in reference^. 

In an incommensurate system all atoms or dimers have 
different geometries and environment. It is therefore im- 
possible to compute all symmetry independent set of t, 
J, e, etc. parameters. Indeed, the atomic positions of 
the present compounds are described as the three dimen- 
sional projection of a periodic four dimensional struc- 
ture (for further details see reference 3 ^). It comes for the 
atomic positions of the two subsystems 

n(p,n x ,n y ,n z ) = r i (p,Q,0,Q)+R l (n Xl n y ,n z )+f (rj(n z )) 

where i refers to the subsystem and p labels the atoms 
in the unit cell, 

n(p, 0, 0, 0) = xo{p)a + y a (p)b + z Q (p)ci 

is the atomic position of the average structure in the 
reference unit cell 

Ri{n x ,n y ,n z ) = n x a + n y b + n z Ci 

is the translation vector of the average structure of each 
subsystem i, and / (rj(p, n z )) is the modulation vector of 
the subsystem i with the periodicity of the other subsys- 
tem : j. 

n(p,n z ) = (z (p) + n z ) — 



and 



f(n(p,n z )) = ^2 A hn (p) sin (2TrnTi(p,n z )) + 

71=1 

Bi, n (p) cos (2wnTi(p, n z )) 

We thus computed the t—J model parameters at 11 dif- 
ferent positions along the ladders. These positions were 
chosen in order to get a good representation of the incom- 
mensurate distortions. The computed values were then 
extrapolated to the whole system using the above crys- 
tallographic description of incommensurate systems. It 
was done using a Fourier's series analysis, as a function of 
the fourth coordinate r, associated with the incommen- 
surate modulations 30 . All results of this work are given 
as a function of the r coordinate of the copper atom (on- 
site parameters) and of the copper atom of lowest x and 
z coordinate for the interaction terms. 



III. INFLUENCE OF THE STRUCTURAL 
MODULATIONS ON THE LADDER SUBSYSTEM 

A. t - J + V model 

Figure |3| displays the hopping and exchange parame- 
ters for the undoped and calcium-doped compounds. One 
first notices that, as was observed on the chains subsys- 

te: 



the variations of the parameters are quite large. 
For instance, on the x — 13.6 ladder legs, the hopping 
modulations reach up to 160 meV. One should however 
point out that the nominal values of both the hopping 
and exchange integrals are much larger for the ladder 
subsystem than for the chain one due to the nearly 180° 
Ou-O-Cu angles. In fact, this is essentially the ligand 
to metal charge transfer integrals of equations Q that are 
sensitive to the structural modulations. In the chain sub- 
system one of them, let us say t a i, is nearly zero (due 
to the nearly 90 c irc angles) and thus very responsive 
to any structural distortions. In the ladder subsystem 
both t a i and tu are large (due to the nearly 180 c irc an- 
gles). The consequence is that the relative variations of 
the parameters are much smaller on the ladder subsys- 
tems. Indeed, even for the hopping on the ladder legs 
of the Chi3.6 compound, the relative variation reaches 
only 22% of the nominal value, while it was 100% for the 
chains nearest-neighbor (NN) hoppings. It can thus be 
expected that despite their large absolute values, hop- 
ping and exchange fluctuations will not be as crucial for 
the ladders electronic structure, as for the chains subsys- 
tems one. Nevertheless one should remember that these 
modulations are not small enough for their effect to be 
negligible. 
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FIG. 3: Effective hopping (a and b) and exchange (c and d) 
parameters as a function of the modulation crystallographic 
coordinate r. Figures a and c correspond to the undoped com- 
pound and figures b and d to the calcium-doped compound 
(x — 13.6). The circles and squares are the computed points 
for the legs and rungs, the solid line are the fitted Fourier 
series. 
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Eccleston et a? 1 reported exchange values extracted 
from neutron scattering experiments, for the undoped 
compound. These values of 130 meV for the ladder legs 
and 72meV for the ladder rungs are to be compared 
with our computed mean values of respectively 143 meV 
and 111 meV. While the order of magnitude are in good 
agreement the ratio between the average rungs and legs 
exchange is larger in our calculations (J±/J\\ = 0.77) 
than in the experimental evaluation (0.55). One should 
however notice that the J±/Ju evaluation from neutron 
scattering does not take into account the four copper 
cyclic exchange parameter. Matsuda et a! 32 found that 
the consideration of this cyclic term in the fit of the neu- 
tron scattering data of the LcieCasOu24:0 '4,1 compound 
notably increases the a J±/J\\ ratio, more in agreement 
with the present results. In addition, the cyclic exchange 
has been evaluated by ab-initio methods for the SrOu20^ 
compound— and found to be 18% of the leg exchange 
value, in qualitative agreement with the neutrons scat- 
tering data of Matsuda et al which evaluate this term to 
15% of the leg exchange. Another argument in favor of 
the quality of our calculations is the consistency between 
the exchange and hopping parameters, computed inde- 
pendently. Indeed, equation ^ yields J±/J\\ — (t±/t\\) , 
while our computed average values gives J±/J\\ = 0.77 

and (Wt||) 2 = 0.75. 

It is interesting to notice that our calculation shows a 
variation of the exchange constant along the ladder legs 
of about 19% of the exchange average value. This result 
compares nicely with the predictions of Schmidt et at^ 
that postulated that the structural modulations in the 
ladders would induce a transfer integral modulation and 
thus an exchange integral one. In addition, they showed 
that such a modulation (~ 16% in their work) of the 
leg exchange values was able to explain the discrepancies 
observed between the Raman spectra of the Srx^Cui^O^x 
and the La§Ca%Cu2AO '41 compounds. 

One can point out that both hopping and exchange 
absolute values increase with the calcium content. This 
is specially true for the ladder rung exchange that goes 
from a mean value of 112 meV for the x = compound 
to 136 meV for the x = 13.6 compound. 

Figure 2] reports the on-site orbital energies, e(r), for 
the x = (c) and x — 13.6 (d) compounds. One sees 
immediately that the variations of these on-site energies 
are large, however much weaker in amplitude than for 
the chain subsystems. Indeed, for the undoped com- 
pound, e varies in a 0.3 eV range while it is 1.2 eV 21 
for the chain subsystem. Similarly for the x = 13.6 com- 
pound the variations range is 1.6 eV for the ladders while 
it is 2.4 eV— for the chains. This results are in agree- 
ment with the hopping and exchange results and are due 
to the facts that (i) the structural modulations are of 
weaker amplitude on the ladder subsystem than on the 
chain one, (ii) the parameters modulation of one sub- 
system, and specifically the on-site orbital energies, are 
essentially due to the structural modulations of the con- 
cerned subsystem itselfii. It is remarkable that these on- 
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FIG. 4: On-site orbital energies differences between NN sites 
(a for x = and b for x — 13.6). The circles correspond to the 
computed values on the legs fragments and the squares to the 
rungs ones. The lines are the fitted Fourier series. Deduced 
on-site energies are presented on figures c for x — and d for 
x = 13.6. The solid and dashed curves correspond to the two 
legs of a ladder. 

site energies modulations are totally correlated with the 
Madelung potential modulation, both for the chain and 
ladder subsystems and both for the undoped and highly 
calcium doped compounds (see figure|SJ|. Indeed, despite 
the fact that the local environments of the ladders and 
chain magnetic centers are very different the computed 
on-site orbital energies exhibit the same dependence to 
the Madelung potential : Ae = 0A2AVMad.Pot. (where 
VMad.Pot. is the Madelung potential on the copper sites 
and delta refers to the NN differences). 




AV MadP0 , (eV) 

FIG. 5: Computed NN-difference of the on-site orbital ener- 
gies as a function of the corresponding Madelung potential 
differences. 

Comparing the e values as a function of the calcium 
doping, one sees that as for the chain subsystem the 
calcium-doped compound presents much larger on site or- 
bital energy fluctuations than the undoped system. This 
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difference is in fact much larger than for the chain sub- 
system since the x = 13.6 versus x — average orbital 
energy ratio is 5.3 for the ladders while it is only 2 for the 
chains. Another important point, that will be discussed 
in more details later on, is the fact that the e variation 
range is weaker than the hopping integral in the x = 
system (0.3 eV versus 0.660 eV), while it is much larger 
in the x = 13.6 compound (1.6 eV versus 0.670 eV) and 
can be expected to induce in this case a hole localization 
on the ladders. 

Another important remark concern the orbital energy 
difference between two sites belonging to the same rung. 
Indeed, for the undoped system this difference is strongly 
modulated, with a NN energy difference larger than be- 
tween NN leg atoms (see figureUK). For the highly doped 
compound however, the two copper sites of a same rung 
are iso-energetic. Indeed, the maximum absolute value of 
the computed orbital energy difference between two sites 
on the same rung is smaller than 1 meV. This energy 
degeneracy witnesses a hidden extra symmetry since the 
two copper sites of a rung respectively correspond to r 
and — t. One retrieves this symmetry on figure^Ji where 
the e(r) function is symmetric around r = 0.5. The ex- 
pected consequence of such a symmetry is the fact that 
the holes present in the ladder subsystem should be de- 
localized between the two sites of the ladder rungs, while 
localized in the c direction. 



IV. CHAIN-LADDER INTERACTIONS 

In this section we will address the possibility of a chain- 
ladder coupling due to the structural modulations. In- 
deed, as mentioned in the introduction, in the highly 
calcium-doped compounds, the chains can be distorted 
in such a way that their oxygen atoms come into an 
apical position of the ladder copper atoms with a small 
Oui a dder-0 c hain distance (2.7 A). In the undoped com- 
pound however, the Cu\ a dder-0 'chain distances always re- 
main larger than 3.1 A. 



tween chain and ladder closest copper sites. One sees 
immediately that unlike what has been supposed in the 
literature, the chain-ladder hopping terms remain weak 
whatever the Oui a dder-0 chain distances. Indeed, both 
the nominal values and the variation ranges (in meV 
units : [—11.7, 1.5] for the x = system and [—18.5,5.5] 
for the x = 13.6 compound) are of the same order of mag- 
nitude for the undoped and calcium-doped compounds. 

We also computed the chain-ladder magnetic couplings 
and the on-site orbital energy differences between the 
chain and ladder systems. The effective exchange were 
always found weaker than 1 meV, that is negligible. 
As far as the chain-ladder orbital energy differences are 
concerned, the results obtained from the present chain- 
ladder fragments are in full agreement with the results 
presented on the ladder subsystems and those from ref- 
erences2ii22 on the chain subsystems. One should point 
out that the present calculations allow us to set the chain- 
ladders relative energy levels, as will be discussed further. 
Let us also notice that the orbital energy differences be- 
tween chain and ladder sites are not proportional to the 
Madelung potential differences on the related sites, un- 
like what was observed for the intra-subsystems orbital 
energy modulations (see figure [SJ. 

One should notice at this point, that counter ions were 
not included in the fragment since their effects were ex- 
pected to be smaller than the contribution due the Ou- 
O-Cu path. Their contributions indeed involve Cki-O- 
Ca/Sr-O-Cu path with Ca/Sr-O distances larger than 
2.4 A. As the obtained values of the interactions are very 
small, their effects could be not totally negligible. Never- 
theless, if the inclusion of the counter ions could slightly 
modify the values of the interactions, it will not change 
the main conclusion which is that the ladder-chain inter- 
actions are largely smaller than the intra-chain or intra- 
ladder interactions. Both hopping and exchange values 
are found more than one order of magnitude weaker. The 
chain and ladder layers can thus be safely considered as 
non-interacting. 
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FIG. 6: Effective hopping between the chain and ladder clos- 
est copper sites as a function of the difference between the two 
copper coordinates along the caxis (Az(Cu)). Circles repre- 
sent the computed points and the dashed lines are a guide for 
the eyes. 

Figure displays the effective hopping integrals be- 



V. DISCUSSION AND CONCLUSION 

One of the most controversial subject about the 
Sri4- x Ca x Ou,2404i family of compounds is the chain to 
ladder hole transfer as a function of the calcium content. 
While it is clearly established that the hole transfer in- 
creases both with the calcium content and the applied 
pressure, the actual amount of transferred holes does not 
reach a consensus, even for the x = system. This sub- 
ject is of importance for the understanding of the system 
properties since (i) it is supposed to be crucial for the 
superconducting phase- and (ii) it is more generally di- 
rectly related to the systems conduction properties since 
the conduction is supposed to occur in the ladder sub- 
system. 

In such a compound, the localization of the electrons 
will mainly arise from (i) the modulation of the on site 
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energies due to structural modulations (ii) the correla- 
tion between magnetic electrons. A full treatment of this 
many body problem in an incommensurate system, is be- 
yond the aim of this work, since it focuses on the role of 
the structural distortions. Moreover, the variations of 
the orbital energies are very large compared to the other 
parameters of this model. In the following, we will thus 
consider only the very first effect of electrons interac- 
tions, by including NN repulsions. This approach will 
already give indications on the role of the distortions on 
the electron localization, and transfer. 

We evaluated the relative energy of the chain and lad- 
der subsystems, responsible for the possible hole transfer, 
from the following parameters : 

• the on-site orbital energies, 

• the NN bi-electronic repulsion in the chain, the NN 
and NNN bi-electronic repulsion in the ladder. 

In references2i*2£, the NN bi-electronic repulsion V was 
computed for the chain subsystem using three centers 
fragment calculations. It was found to be almost inde- 
pendent of the structural modulations, temperature and 
calcium content. This can be understood by the fact that 
the NN repulsions are essentially dependent on the Cu- 
Cu distances and not on the Cu-O-Cu angles that dom- 
inate the modulations of the other effective interactions. 
We thus evaluated the NN repulsions on the ladders sub- 
system using the average computed chain value of 0.6 eV 
and the standard Ohno formula^ for the distance depen- 
dence 



V(R) 



Vn 



1 + R/a 



where R is the inter-site distance, do is the Bohr radius 
and Vo a R = effective constant. Extracting Vo from 
the chain repulsions we get Vo — 3.7 eV. We thus find 



for the intra-ladder repulsions : Vi 
0.44 eV, and V d 
on a plaquette), 



< !l 



v, 



rung 



0.32 eV (diagonal repulsion 



for the inter-ladder repulsions : V 



eg-leg 



0.6 eV. 



Using these approximations and the computed orbital en- 
ergies we determined the relative energies of the chains 
and ladders as a function of t (see figure 0. Since the 
ladder hole-doping remains small we considered that the 
repulsions where acting on all ladder sites. As a matter 
of comparison we computed the mean (periodic approx- 
imation) relative energy between the ladder and chain 
sites using the same method. We found for the x — 
undoped compound the chains about 0.5eV above the 
ladders while for the x — 13.6 compound the chain and 
ladders subsystems are quasi-degenerate. For the real 
compounds, one can see on figure that while for the 
x = compound the ladder energy curve is always about 
0.65 eV lower than the first empty sites of the chains, for 
the x = 13.6 compound, the most energetic ladder sites 
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FIG. 7: Chain and ladder energies as a function of Trains- 
The site energies were evaluated as the sum of the orbital en- 
ergies and the repulsion terms with neighboring atoms when 
the latter have lower orbital energies and are thus occupied 
prior to the considered one. The ladder fourth coordinate 
ladder has been rescaled to 2c chain /c ladder , so that the rela- 
tive r variation range of the two subsystems is proportional 
to the number of sites per f.u. in each sub-system. The hori- 
zontal line represents the Fermi level. 



are at the Fermi level. These results are coherent with 
a complete hole localization on the chain sub-system for 
the undoped compound and a small hole transfer to the 
ladders for highly doped systems. It seems however diffi- 
cult to extract a precise value of the chain to ladder hole 
transfer (for the x = 13.6 compound). Indeed the hole 
transfer can be expected to be sensitive to very small 
energetic variations, since the Fermi level is located in a 
region with large density of states, both for the chain and 
ladder subsystems. For the ladders subsystem, Ep is lo- 
cated on an energetic plateau. For the chain subsystem, 
Ep crosses the energy curve, unlike what happens in the 
undoped system where the Fermi level is located in an 
energetic gap. This may explain the large experimental 
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range of values obtained for this hole transfer. 

Another consequence of the Fermi level localization on 
a ladders energetic plateau, is that the holes are local- 
ized on (quasi) iso-energetic rungs — while they are lo- 
calized on very energetically different sites in the chain 
subsystem. These rungs are not first neighbors, but in 
most cases second and in some cases third neighbors. The 
holes can thus be considered as evolving on a subset of the 
ladders with (quasi) iso-energetic sites. In such a model, 
the rungs belonging to the real ladder located in between 
two hole-supporting rungs can be considered as bridging 
ligands. The resulting effective hopping and exchange in- 
tegrals between the hole-supporting rungs are thus mod- 
ulated according to the number and the characteristics of 
these in-between rungs. In such a model, these effective 
interactions between hole-supporting rungs can be ex- 
pected to be rescaled to much lower values than the NN 
ones on the real ladder, and thus much lower values than 
the on-rung ones (\tf® g \ < |t rung | and |Jf*| < lining I)- 

We would like to conclude this work on a few opened 
questions on the super-conducting phase. First, the on- 
site energy variations are quite-large on the calcium- 
doped system, both for the chains and ladders subsys- 
tems. These modulations strongly increase as a function 
of the chemical pressure. Whether the applied pressure 
would act along the same line, it would induce an even 
larger distortion of the chains and ladders subsystems. 
Such a phenomenon would be quite incompatible with 
the usual super-conductivity theory in the ladder systems 
since it would induce an even larger modulation of the 
ladder on-site orbital energies. Another discrepancy be- 
tween the present results and the super-conductivity pre- 
dictions is that the super-conductive pairing is supposed 
to occur for (i) low ladder hole doping, (ii) 2\J rung \ > 
\Jieg\ (both are verified in the present case) and (hi) 
l^rimgl > l^rwngl which is far from being verified since 
1^™^^! > 4| Jrungl- One thus wonder either what type of 
structural distortions under applied pressure could cor- 
rect these features hindering the super-conductivity or 
whether the observed super-conductivity arise from an- 
other origin that the usually accepted homogeneous lad- 
der t — J model. 
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Appendix 



The computed results were fitted as a function of the 
fourth crystallographic coordinate r, using a Fourier se- 
ries, according to the following expression 









a n cos [2 


im(r — ip 


n)] 


(2) 








n 


















(a) 








x = 




El 


A£ r 


ti 


t r 


Jl 




do 




0.0 


0.0 


-659.2 


-572.2 - 


143.9 


-111.6 


ai 




153.0 


204.3 


28.3 


12.4 


-13.7 


2.5 


ai 




8.9 


3.4 


2.0 


-5.1 


-1.0 


-2.5 


Vi 




0.819 


0.875 


0.655 


0.750 


0.904 


0.750 






0.786 


0.688 


0.089 


0.500 


0.982 


0.500 










(b) 








x = 13. 


6 


El 


Ae r 


ti 


t r 


Ji 




a 


0.0 


0.0 


-672.5 


-624.7 


-159.3 


-136.4 


CL\ 


-811.4 




73.9 


11.5 


-0.2 


8.0 


02 


-209.5 




-11.1 


4.0 


-10.5 


4.5 


0,3 










4.4 




<fl 


0.500 


0.000 


0.284 


0.001 


0.284 


0.995 




0.004 


0.000 


0.529 


1.000 


0.537 


0.997 












0.972 




TABLE I: Analytic fit of the t - 


J model, 


a) for the x = 



undoped compound, b) for the x = 13.6 calcium-doped 
compound. All energies are given in meV. 



Only terms with a non negligible contribution to the se- 
ries were retained. 

The results are summarized in tables [IJi for the x = 
and [J) for x = 13.6. 
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